function dni = DNI(D,t)
%法向直接辐射辐照度计算
%参数 %D为距离春分日的天数 t为当地时间（真太阳时）24小时制

global H;%海拔
H=3; %3000m

global lon %(东)经
global lat %(北)纬
lon=98.5;
lat=39.4;
lat=deg2rad(lat);
sin_sigema=sin(2*pi*D/365)*sin(2*pi*23.45/360);
sigema=asin(sin_sigema);%sigema为太阳赤纬角,其中asin是求反正弦弧度值
omega=pi/12*(t-12);%太阳时角计算


sin_alpha_s=cos(sigema)*cos(lat)*cos(omega)+sin(sigema)*sin(lat);

%DNI计算公式系数
a=0.4237-0.00821*(6-H)^2;
b=0.5055+0.00595*(6.5-H)^2;
c=0.2711+0.01858*(2.5-H)^2;

G=1.366;%太阳常数

dni=G*(a+b*exp(-c/sin_alpha_s));

end

